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USE  OF  INCOHERENCE  TO  PRODUCE  SMOOTH 
AND  CONTROLLABLE  IRRADIATION  PROFILES 
WITH  KrF  FUSION  LASERS 

INTRODUCTION 

Direci-drive  laser  fusion  requires  high  intensity  laser  light  that  uniformly  illuminates  a  spherical 
target  surface.  Theoretical  studies  have  shown  [1-3]  that  acceptable  spherical  uniformity  can  be 
achieved  with  direct  illumination  by  overlapping  a  limited  number  (>20)  of  focused  beams,  provided 
that  each  individual  beam  profile  is  smooth  and  reproducible.  Earlier  efforts  to  obtain  such  profiles, 
however,  had  been  frustrated  by  the  inherent  imperfections  in  high  power  multistage  lasers.  The 
cumulative  effect  of  small  amplitude  and  phase  aberrations  introduced  by  each  optical  element  of  a 
multistage  laser  is  to  produce  large  random  aberrations  in  the  output  beam.  In  the  usual  configuration, 
where  the  target  is  placed  in  the  quasi  near-field  of  the  focusing  lens,  these  aberrations  tend  to  produce 
large  random  intensity  nonuniformities  at  the  target  surface. 

One  possible  solution  to  this  problem  is  a  technique  called  Induced  Spatial  Incoherence  (ISI) 
[4.51.  In  this  scheme,  spatially  incoherent  light  is  created  by  propagating  a  laser  beam  of  broad  spectral 
bandwidth  through  a  pair  of  echelon  structures  that  impose  a  different  time  delay  at  each  step.  If 
the  delay  increments  are  chosen  somewhat  larger  than  the  optical  coherence  time  tr  —  1/Av,  the 
beam  is  sliced  into  an  array  of  mutuaily-incoherent  square  beamlets.  Each  of  these  will  independently 
focus  to  the  same  diffraction  profile  [the  sinc2(jr)  sine2 O')  function]  of  width  f\/d,  provided  that  its 
initial  width  d  is  small  in  comparison  to  the  transverse  scalelength  s„  of  the  incident  beam  aberration. 
(Here,  /  is  the  focal  length  of  the  lens,  and  X  is  the  mean  optical  wavelength.)  The  transient  interfer¬ 
ence  pattern  produced  by  superposition  of  these  beamlets  will  evolve  randomly  in  times  of  order  tc. 
The  target  will  therefore  ignore  this  rapidly  shifting  structure,  and  respond  only  to  the  time-averaged 
intensity  <I(x)>  <*  sinc20c)  sine2 O')  if  its  hydrodynamic  response  time  th  satisfies  th  >  >  *r.  For 
example,  an  optical  bandwidth  Av  —  30  cm-1  (easily  achieved  in  Ndiglass  or  KrF  lasers)  provides  rc  = 
l  psec.  whereas  th  is  typically  =  l  nsec  in  large  high-gain  pellet  designs. 

In  its  present  application,  the  ISI  technique  requires  the  echelons  to  be  placed  at  the  output  of  the 
laser.  This  restriction  stems  from  self-focusing  effects  in  glass  lasers,  and  from  the  necessity  of  main¬ 
taining  spatial  coherence  in  any  harmonic  conversion  crystals.  The  near-field  nonuniformities  associ¬ 
ated  with  an  array  of  beamlets  would  seed  self-focusing  if  one  attempted  to  amplify  those  beamlets  in  a 
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multistage  glass  laser.  [6]  where  the  nonlinear  phase  shifts  are  typically  B  >  5  radians.  By  placing  the 
echelons  at  the  output,  one  can  control  seif-focusing,  and  thus  maintain  an  acceptable  degree  of 
transverse  beam  uniformity  over  distance  —  d.  This  configuration  wouid  require  a  large  number  of  op¬ 
tically  coated  echelon  steps  operating  at  high  optical  fluence  levels  if  it  were  used  in  a  fusion  reactor 
with  large  apertures  and  multiple  beam  lines.  For  example,  a  recent  conceptual  design  study  for  the 
Sirius-M  reactor  [71  concluded  that  in  order  to  use  ISI  at  X  *  1/4  ^m,  one  would  require  240  steps  in 
each  transverse  direction  for  each  of  32  drive  beams.  Another  issue  raised  by  this  ISI  configuration  is 
that  of  efficiency.  Approximately  15%  of  the  energy  at  the  focal  plane  will  diffract  into  sidelobes,  and 
most  of  that  energy  would  have  to  be  discarded  in  order  to  achieve  good  illumination  uniformity. 

Although  self-focusing  remains  a  serious  problem  in  glass  lasers,  it  is  far  less  important  in  KrF 
systems  [8-131.  where  intensities  are  typically  <  10  MW/cm2,  and  the  amplifying  medium  is  gaseous. 
A  reduction  of  the  nonlinear  phase  shift  to  <  1/2  radian  could  eliminate  the  necessity  of  placing  the 
echelons  at  the  laser  output,  thereby  opening  the  way  for  several  possible  improvements.  For  example, 
instead  of  the  reactor  configuration  discussed  above,  one  could  produce  the  beamlets  by  a  single  pair  of 
echelons  at  a  low  power  stage,  spatially-filter  them  to  eliminate  the  sidelobes,  then  optically-relay  them 
through  each  of  the  main  amplifier  chains.  As  long  as  these  beamlets  remain  small  in  comparison  to 
the  transverse  scalelength  of  the  aberrations  (e.g.,  due  to  passive  optics,  turbulence,  and  nonuniformi¬ 
ties  in  the  amplifier  excitation),  they  can  focus  to  the  sine2  profile  without  any  sidelobes  at  the  target. 
One  issue  raised  by  this  scheme  is  the  fact  that  the  spatially-filtered  beamlets  would  have  Gaussian-like 
intensity  profiles  that  remain  stationary  within  the  amplifier  apertures;  hence,  there  is  a  possibility  that 
they  (and  therefore  the  target  beam  profile)  could  be  affected  by  gain  saturation. 

This  paper  describes  an  alternative  ISI  scheme  shown  in  Fig.  1,  in  which  the  echelons  are  com¬ 
pletely  eliminated.  The  concept,  which  will  be  referred  to  as  echelon-free  ISI,  is  basically  an  image  pro¬ 
jection  technique  that  projects  the  desired  spatial  profile  onto  the  target  via  the  laser  system,  using 
partially-coherent  light.  The  information  required  to  reproduce  this  profile  is  transported  through  the 
laser  by  a  multitude  of  small  coherence  zones,  rather  than  any  large  whole-beam  structure.  Thus,  the 
coherence  zones  play  essentially  the  same  role  as  the  beamlets  do  in  conventional  ISI,  except  that  their 
near-field  intensity  profiles  are  not  "frozen  in"  by  any  echelon  steps.  A  large,  spatially-incoherent  beam 
of  amplitude  £0(x,r)  is  generated  by  broadband  amplified  spontaneous  emission  (ASE)  in  a  mirrorless 
or  spatially  multimode  laser  oscillator.  Its  central  portion  illuminates  a  transparency  or  variable-density 
absorber  VDA.  whose  spatial  transmission  function  F(x)  (of  effective  aperture  D0)  is  the  smoothly- 
varying  irradiation  profile  required  by  the  target.  [2.3]  By  choosing  the  oscillator  aperture  sufficiently 
large,  and  placing  it  close  to  VDA,  one  can  satisfy  the  following  conditions:  (i)  The  light  within  D0 
must  remain  incoherent,  with  a  broad,  smooth  angular  power  spectrum  y^B)  of  width  Ad>>\/Z?o. 


DESIRED  IRRADIATION  PROFILE  F(r)  IS  IMAGED  ONTO 
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(ii)  To  a  good  approximation,  it  must  be  statistically  homogeneous;  i.e.,  both  yo(0)  and  the  time- 
averaged  incident  intensity  </0>  -  < !  £0(x)l2  >  must  remain  nearly  uniform  across  D0.  Any  devia¬ 
tion  from  condition  (ii)  is  likely  to  arise  from  the  largest  angie  modes  of  yo(0),  which  may  experience 
some  gain  nonuniformity  within  the  ASE  oscillator.  These  modes  can  be  readily  eliminated,  however, 
by  slightly  overfilling  the  aperture  Dx  at  the  iaser-chain  entrance  pupil  EP,  as  shown  in  Fig.  1.  If  condi¬ 
tions  (i)  and  (ii)  are  satisfied,  then  the  beam  £0'(x,r)  transmitted  through  VDA  can  be  described  by 
the  spatial  autocorrelation  function 

r0(I,7)  =  <£o(jc)£o*Cc')>  =  </0>f(l)y0(?).  (1) 

Here.  X  =  y(3?  +  x')  and  7  =  x  -  x '  define,  respectively,  the  mean  position  and  separation  of  the 

transverse  coordinates,  while  yo(f)  (which  is  the  Fourier  transform  of  yo(9)\  is  nonvanishing  only 
within  a  coherence  zone  d$  >  \Q<  <  Z)q  (E.g.,  see  Fig.  1).  The  correction  term  to  the  approxima¬ 

tion  in  (1)  is  of  order  dll Dq.  which  is  negligible  for  the  ratios  d0  <  Do/40  envisioned  for  actual  KrF 
systems. 

Because  VDA  and  EP  are  located  at  opposite  focal  planes  of  lens  L\,  the  amplitude  £f  (x.  t)  at  EP 
will  be  the  Fourier  transform  of  £0(x,  t).  The  autocorrelation  function  <  £t  (jc)  £:  *(F)>  at  EP  then 
becomes 

r,(.r.7)  -  </|(j?)>  £(7).  (2) 

where  <Ix(X)>  -  <|£,(7)l2>  «  y0(7//t)  is  the  average  intensity,  which  may  vary  slowly  across 
D\,  and 

£( 7)  -  f  F(X')  exp  (-  i  X'  •  7)  d1X\  F{ 0)  -  1 .  (3a) 

—  0  for  (71  >  dx  =  K/i/D^  <<  Dx  (3b) 

is  the  Fourier  transform  of  the  absorber  function,  whose  coherence  width  is  d\  —  \f\l  Do  «  D\.  For 
example,  a  Gaussian  profile 

F(X)  -  (1/ir  Rl)  exp  (-  l.?l2/£02)  (4) 

of  Me  diameter  0O  -  2i?o  would  transform  to  £(f)  -  exp  (-  if  l2/rf),  with  a  coherence  zone  of  diam¬ 
eter  d\  -  2r,  -  (4/ir)  x/i/Dq.  Expressions  (3a.  b)  show  clearly  that  the  intensity  and  phase  informa¬ 
tion  needed  to  reproduce  F(X)  is  now  contained  in  the  small  coherence  zones  of  width  </,.  rather  than 
any  large-scale  beam  structure  of  width  ~DX  ^Fig.  1). 


Beyond  EP.  the  light  is  relayed  through  the  aberrated  multistage  amplifier  chain  to  produce  an 
output  amplitude  E2(x,  t)  at  the  focusing  lens  L  2.  (The  effects  of  aberration  are  illustrated  in  Fig.  1  by 
the  non-symmetric  appearance  of  the  average  intensity  <  i£i(T)  l2> .)  The  average  intensity 
<1  £r  (x)  l2>  at  the  target  is  then  given  by  the  general  expression 

</r  < *)>  -  ~T  Sdlx'S  d1*'  r2  <*'•«'>  exP  "  i  4-  *  ■  T  .  (5) 

Si 

where  T :  UT,£  )  is  the  autocorrelation  function  at  L,.  In  the  ideal  case,  where  self-focusing  is  negligi¬ 
ble,  the  coherence  zones  at  L2  will  convey  essentially  the  same  information  as  those  at  EP,  provided 
that  they  remained  small  in  comparison  to  the  local  apertures  and  aberration  scalelengths  s„,  as  they 
propagated  through  the  laser  system.  The  only  change  will  be  the  magnification  factor 
M  -  D-J  D{  —  d-J  d\  due  to  the  expanding  apertures.  Thus,  the  output  beam  will  focus  to  the  desired 
average  intensity  profile  of  width  DT  =  \fjd2  -  Z)0 

< 7r (jc ) >  -  (A//i//:)2  <P>  £(-  ) .  (6) 

where  <P>  is  the  average  output  power.  The  transient  interference  structure  (of  typical  size 
dT  -  \ f  -J Di  >  d<rf\JfxM)  will  average  out  over  many  coherence  times,  just  as  it  did  in  the  original 
ISI  technique.  In  the  remainder  of  this  paper,  we  will  examine  the  perturbing  effects  of  linear  and  non¬ 
linear  aberration,  and  present  preliminary  experimental  results  demonstrating  that  the  technique  can 
produce  smooth  target  beam  profiles  even  under  non-ideal  conditions. 

Because  this  scheme  uses  image  projection,  it  bears  some  resemblance  to  techniques  that  use  exci- 
mer  lasers  for  projection  photolithography  (141.  The  essential  difference,  however,  is  that  the  profile 
Fix)  lacks  the  fine  structure  normally  demanded  in  high  resolution  lithography;  consequently,  it  does 
not  require  large  coherence  zones  that  would  be  susceptible  to  aberration  in  the  laser  system.  By  con¬ 
trast.  projection  photolithography  requires  considerable  care  to  minimize  aberration  in  the  lenses  and 
image  projection  path. 

Echelon-free  ISI  is  compatible  with  angular  multiplexing  techniques  proposed  for  excimer  fusion 
drivers  (9-131.  For  example,  a  5  mm  coherence  zone  (appropriate  for  a  30  cm  aperture)  would  have  a 
50  jitad  diffraction  angle  at  X  -  0.25  jam,  whereas  the  angular  separation  of  the  multiplexed  beams 
would  typically  be  >1  mrad.  The  beam  at  EP  need  not  be  exactly  imaged  onto  lens  L->:  in  fact  a  sig¬ 
nificant  amount  of  free  propagation  would  be  tolerable  within  the  amplifier  system.  The  only  require¬ 
ment  is  that  any  spreading  of  the  coherence  zones  due  to  diffraction  must  also  remain  small  in  com¬ 
parison  to  the  apertures  and  aberration  scalelengths.  Thus,  the  technique  would  not  require  image 
relaying  at  the  large-aperture  final  amplifier,  or  even  in  the  temporal  decoding  paths  [9-131  that  tollow 
it.  For  example,  the  diffraction  length  :d  =  <i2/\  of  a  5  mm  coherence  zone  *ouid  be  ~100  meters 


(300  as).  Image  relaying  would  be  required  in  most  of  the  amplifier  chain,  including  the  temporal 
encoding  paths  near  the  input.  For  the  small  apertures  near  the  input,  however,  a  separate  telescope  in 
each  encoding  path  would  be  a  minor  price  to  pay  for  controllable  high  quality  beams  on  target. 


ABERRATION  EFFECTS 

Using  a  simplified  model  of  multiplexed  excimer  lasers,  we  have  carried  out  a  detailed  perturba¬ 
tion  analysis  to  evaluate  the  effects  of  linear  aberration,  self-focusing,  amplifier  saturation,  and  diffrac¬ 
tion  on  the  focal  profile  with  echelon-free  ISI.  This  section  will  outline  the  assumptions  and  results  of 
that  analysis;  a  more  detailed  derivation  will  be  given  in  a  later  publication.  A  perturbation  treatment 
should  be  appropriate  for  this  problem  because  the  effects  due  to  aberration  must  remain  small  if  the 
technique  is  to  be  at  all  useful.  In  treating  the  individual  aberration  mechanisms,  the  key  simplifying 
assumption  is  the  thin  window  model,  which  ignores  the  effects  of  diffraction  in  the  locality  at  which 
the  aberration  is  created.  This  can  be  justified  by  the  fact  that  most  aberration  occurs  in  optical  ele¬ 
ments  short  in  comparison  to  the  diffraction  lengths  zd  -  d2/\  associated  with  the  coherence  zones, 
(e.g..  recall  the  :d  estimate  given  above.) 

Linear  aberration  can  arise  from  turbulence  and  nonuniform  excitation,  or  from  passive  optical 
elements  such  as  lenses,  mirrors,  and  apertures.  In  the  thin  window  approximation,  it  is  described  by  a 
nonuniform  intensity  transmission  T(x)  and  phase  shift  <*fjc) ,  which  relate  the  incident  and-aberrated 
amplitudes  by  the  expression  E'(x,r)  -  £(3r.r){r(x))1,2exp[/<6(x)].  The  corresponding  correlation 
functions  therefore  satisfy 


r(f,?>  -  ru,*)  (m)  ronil/2  exP  u<t>m  - 


where  x  -  X  -!-  yf  and  X  -  yf.  Because  If)  is  limited  to  coherence  zones  <i  small  in  com¬ 
parison  to  the  aberration  scalelength  s3,  one  can  approximate  this  expression  by  a  second  order  expan¬ 
sion  in  the  small  quantity  if  \/sa-~d/sa\  i.e., 

r*  (.?.?)  =  nxj)  nx)  u  +  /?  •  vd>(x)  -  If  •  ( (v«)  (vt)  -  ^  vvu„t)  l  •?),  (i) 
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where  V  =  d/dX  operates  in  the  transverse  plane. 


Self-focusing  is  a  potentially  serious  limitation  because  it  can  affect  the  beam  on  scaielengths  com¬ 
parable  to  the  coherence  zones.  The  main  problem  arises  from  those  contributions,  such  as  the  third- 
order  electronic  susceptibility.  [6!  whose  nonlinear  refractive  index  can  respond  "instantaneously”  to 


•.*  -. 

■  •r  s.  <  --  -• 


the  stochastic  intensity  fluctuations  (i.e.,  on  a  time  scale  <  r,).  In  the  thin  window  model,  this  contri¬ 


bution  will  generate  an  instantaneous  nonlinear  phase  shift  rji£(x.:)  I2,  where  tj  -  (8ir//»oc)y  n2kL, 


/to  is  the  linear  refractive  index,  and  L  is  the  thickness  of  the  medium.  The  transmission  through  this 


medium  is  then  described  by  the  relation 


r  (Jr,?)  -  <E(x,t)Em(x',t)  exp  [/7)|£(jf,r)l2  -  /t,|£(x',/)!21>. 


where  (X,£)  are  related  to  (x,x')  in  the  usual  way.  Because  the  incident  tight  is  chaotic,  the  ampli¬ 
tudes  £(x,c)  and  E(x\t)  should  satisfy  Gaussian  statistics;  it  is  therefore  permissible  to  expand  the 
exponential  and  apply  Gaussian  factorizations  to  the  resulting  amplitude  products.  Carrying  this  out  to 


second  order,  and  recalling  that  f  is  a  small  parameter,  one  obtains  the  approximate  result 


rur,7)  a.  r(J.?)  {1  +  2/*  ■  v<b(X)>  -  211  -  lyur,£)l2i  <£un>2}. 


where  y(X,£)  =  f  (X,(-)/ <  I(X)>  ,  and 


<B(X)>  =  tj< I £( J) |2>  -  t,  </(X)> 


is  the  average  nonlinear  phase  shift,  which  must  be  kept  as  small  as  possible  (  <B>2  <  <  1).  For 
fused  quartz  (no  =»  1.5,  rt2  =  10_u  esu  [61),  Eq.  (9)  gives  <B>  =  7  x  10~5  <l(MWIcm1)> 

■  L  (cm)  at  A  -  0.25#im. 


The  aberrating  effects  of  gain  saturation  depend  upon  how  rapidly  the  stochastic  intensity  fluctua¬ 
tions  I(x,t)  —  <I(x)>  occur  within  the  relaxation  time  r  of  the  excimer  population.  For  fluctuation 
times  -  i;<  <  t  in  a  homogeneously  broadened  laser  medium  (e.g.,  KrF).  [91  the  large-signal  gain 
G(x)  will  respond  only  to  the  average  intensity  </(*)>.  The  relationship  between  the  incident  and 
amplified  correlation  functions  thus  follows  from  the  elementary  expression  E'(x,t)  —  lG(x)]in 
E(x.t).  Because  G(x)  —  G(<l(x)>)  varies  slowly  within  the  coherence  zones,  it  can  be  expanded  in 
the  same  manner  as  the  linear  transmission  T{x)  treated  previously;  i.e.. 


rd,£)  =  H*,?)  G(X)  l  1  +  -£-?•  (V  V  In  G)  -71. 

o 


Using  a  perturbation  treatment  of  the  laser  amplifier  equations,  one  can  show  that  the  correction 
term  due  to  intensity  fluctuations  with  nonvanishing  values  of  tc/r  is  approximately 


a  r ■(*.?)=»  -  (*,m  -  g2(X)  I  y(x,o  I2 1  r  (.r,,c)  tc/r. 


where  yCV.f)  is  the  reduced  correlation  function  defined  previously,  and  g  ( ( A' ) .  £2<,Y)  are  slowiv- 
varying  functions  satisfying  i  £,  |  —  !g2l  —  1-  If  this  term  were  large,  it  could  be  an  additional  source 
of  small-scale  aberration  due  to  i  y(X.O  I2;  however,  for  typical  values  of  4  —  1  ps  lAv  -  30cm”1) 
and  t  —  2  ns,  it  remains  negligible. 

Eximer  Laser  Model 

The  aberration  terms  can  now  be  incorporated  into  a  simple  laser  model  to  provide  useful  expres¬ 
sions  for  the  target  beam  perturbation.  Aberration  due  to  the  amplifier  chain  alone  is  modeled  by  a 
single  thin  window  Wx  of  nonuniform  gain  G(x)  and  phase  <f>x(.x)  imaged  at  or  near  the  entrance  pupil 
EP.  Thus,  the  autocorrelation  function  r'((Af,f)  at  the  output  of  the  laser  chain  can  be  related  to  the 
appropriately-scaled  input  function  T  x(X/My  £/M)/M2  [Eq.  (2)]  by  expression  (7),  with  the  replace¬ 
ment  T{X)  —G(X)  and  <6(T)  —  <6i(A0.  [In  general,  the  nonuniformity  of  G(X)  may  include  the 
effect  of  apodization  in  addition  to  those  of  nonuniform  excitation  and  losses  within  the  amplifiers.] 
The  single-window  approximation  is  justified  by  the  fact  that  most  aberration  is  likely  to  occur  in  either 
the  amplifier  modules,  which  are  normally  imaged  onto  EP,  [11-13]  or  in  optical  elements  (e.g..  tele¬ 
scope  lenses)  that  lie  in  the  quasi  near-field  of  those  modules.  Self  focusing  effects  are  expected  to 
remain  negligible  within  the  amplifier  chain  because  of  the  low  intensities  required  in  excimer  lasers. 
For  example,  at  typical  KrF  output  intensities  <15  MW/cm2,  [9]  expression  (9)  gives  <B>  <10~3 
L(cm)  in  fused  quartz:  at  these  intensity  levels,  even  a  total  nonlinear  pathiength  of  50  cm  would  gen¬ 
erate  only  <B>  <0.05  radian. 

A  second  thin  window  W2y  which  would  be  located  near  lens  L2,  models  linear  and  nonlinear 
aberration  in  the  focusing  optics  and  delay  path  (required  for  temporal  decoding  [9-13])  beyond  the 
amplifier  chain.  Diffraction  in  this  delay  path  is  modelled  by  allowing  a  free  propagation  distance  :x2 
between  the  windows;  the  correlation  function  IMAf.^)  for  light  incident  at  W2  can  then  be  related  to 
r'1(A'.$)  by  paraxial  diffraction  theory  [15]: 

r;(£,7)  -  1  f d2x-  f d2?  r,<r,f)  exp  [/  —  (r  -  x)  ■  (?'  -  ?)]  (in 

2Jr-12  J  -12 

In  the  spirit  of  the  perturbation  treatment,  this  integral  is  evaluated  only  up  to  second  order  in  the 
small  quantity  dx2Jsa,  where  dx2  =  X  zX2)d2  represents  the  diffractive  blurring  of  a  coherence  zone  of 
initial  width  d2  -  Mdx.  The  focusing  optics  and  delay  paths  are  more  likely  to  introduce  phase  aberra¬ 
tion  <6:(x)  (e.g.,  due  to  thermal  gradients),  rather  than  any  significant  transmission  nonuniformity; 
hence,  the  transmission  T  <  1  will  be  treated  as  uniform  in  I*';.  Combining  this  consideration  with 
the  nonlinear  perturbation  terms  derived  previously,  one  obtains  at  the  output  of  W2, 


r:(.r.£)=  r2u,f)  t  U  +  /  f  •  va,  m  <v«*2)  (v<62>  • « 

+  2/?  ■  V<3(X)>  -  2  [1  -  |y2<X?)l:]  <BIX)>2],  (12) 

where  y2(A\£)  —  r:  (X,§)/<  /2(.V)  >,  and  <B(X)>  is  given  by  Eq.  (9)  with  <HX)>  — *  <  I  AX)  > . 
The  average  target  beam  intensity  <IT(.x)>  can  then  be  evaluated  by  substituting  f2(.Y,  f)  into  Eq. 


Combining  the  above  considerations,  and  retaining  terms  only  up  to  second  order  in  the  small 
quantities  d-,j s3,  d12/sa  and  <B>,  one  can  finally  evaluate  the  perturbed  target  beam  profile 
<IT  (x)>  from  Eq.  (5).  In  terms  of  the  average  output  power 


r  -  <Il(X/M)>  , 

:p>  -  J  T  GW)  - 1 - 5 - d2X' 


and  the  weighted  transverse  spatial  average 


\Q\av  =  f  T  GW)  <h  Q(X')  d2X\ 


which  is  defined  for  any  well-behaved  quantity  Q(X'),  one  obtains 

<IT  (x)>  -  </f0)  (3c)  >  +  <i;u  (3c)  >  +  </f2>  (x) >  +  </f‘vt>  (3c)  > 


where 


c/f0’  (3c)  >  a  (A//,//2)2  <P>  F  (-  Xlf\ 


</fu  C x)>  =-(/j/fc)  {V  •  V  </^0,(3c)>. 


<Ii»  (3c)>=  y  WJk)2  ((Vm'<6)  (V, <6) 


~  jV/V/Zb  (TO  7„  </fD)  (3c)  > 


-  l:\Jk)  [V/V/«A2  U  U„  <//u)  (3c)>] 


</fVL)  <x)>  =  -  2 (f,Jk)  (V'  <B>  \AV  ■  V  </r  ( x)> 
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-  2(  <a>2  ).<„  [  </f01  (x)>  -  (A//,//;)2  <P>  tf(x)  ]. 


ff(x)  =  (-^-)2  f  I  F  (?)  i2  /•  (f)  exp  (-#-7-  .Wx  •  f)  </2?, 


and  <f> ( .V ')  =  (^1  (X')  +  <£2  (AT)  is  the  total  linear  phase  aberration.  [In  Eq.  (16c),  the  repeated 
indices  p,  o  imply  a  summation  over  the  two  transverse  dimensions.]  The  zeroth  order  term  (16a)  is 
identical  to  Eq.  (6).  Expressions  (16b)  and  (16c)  are,  respectively,  the  first  and  second  order  perturba¬ 
tions  due  to  linear  aberration,  and  (16d),  which  is  entirely  in  second  order,  describes  nonlinear  effects. 


In  (16d),  the  first  term  accounts  for  beam  steering  due  to  a  combination  of  nonlinear  refraction  and 
nonuniform  intensity,  while  the  <B> 2  term  describes  self  focusing  in  the  vicinity  of  L2,  which  results 
in  profile  broadening  at  the  target. 


The  small  expansion  parameters  may  be  seen  more  clearly  in  these  expressions  by  noting  that  the 
operators  V  '  =  3/0 A" and  V  =  3/0x  can  introduce  factors  of  order  l/sa  and  l/Dr,  respectively,  where 
DT  -  D0  f-J f\M  is  the  target  beam  diameter.  Recalling  the  relationship  DT  -  Kfi/d2  between  Dr  and 
the  coherence  zone  diameter  d2,  one  can  see,  for  example,  that  expression  (16c)  will  contain  the  fac¬ 
tors  di/$a  and  d2  d{2/s2,  while  (16d)  contains  the  factors  <B>  dij sa  and  <B>2.  Contributions  of 
order  d?2/s2  and  <B>dl2/sa  vanish  under  the  conditions  assumed  here. 


It  is  instructive  to  relate  the  linear  aberration  coefficients  to  the  angular  divergence  and  focal  spot 
parameters  that  the  laser  would  produce  (in  the  absence  of  nonlinear  aberration)  with  a  plane  wave 
source  at  EP.  In  the  geometrical  optics  limit,  this  relationship  can  be  seen  immediately  by  noting  that 
?(,?')  =  (1  Ik)  V '  <b(X')  is  the  refraction  angle  that  an  incident  axial  ray  would  acquire  as  it  pro¬ 
pagates  through  the  aberrated  system.  The  centroid  of  the  aberrated  (but  spatially  coherent)  focal  dis¬ 
tribution  would  then  be  given  by  the  intensity-weighted  average 


=  “  (/;/fc){V  <b(X  )].|i'. 


which  is  the  coefficient  in  </ru(x)>.  If  the  laser  could  be  precisely  aligned  onto  the  target,  then 
[xc)Ay  would  vanish  identically;  in  practice,  however,  the  random  angular  "jitter"  80(A  )  due  to  micro- 
turbulence  8<6(x ')  in  the  amplifiers  and  propagation  paths  will  quickly  lead  to  a  small  residual  misalign¬ 
ment  [xe\4y  /:!80)^k  (Gradual  misalignment  due  to  large-scaie  thermal  gradients  would  presumably 
be  corrected  by  an  automatic  alignment  system.)  For  the  second  moment  \xc  xc}Ay,  geometrical  optics 
would  give  f\  10^1  tv.  however,  a  more  complete  treatment  including  diffraction  yields  the  result 


V 


Sj 

! 


{xc  X'}Ay  -  (/2/*)2{  (V'<6)  (V'<6)  -  jV'V'/rt  ( TG)  }iV.  (18) 

which  accounts  for  the  first  two  terms  of  </f2,(x)>.  The  final  term  in  </f2)  ( x)>  (which  is  the  only 
one  that  involves  the  free  propagation  path)  has  the  coefficient  (zyi/k)  {V '  V ' <62j.,^.  and  arises  from 
the  combined  effects  of  diffraction  and  phase  aberration  beyond  the  laser  chain.  This  expression 
represents  a  phase  curvature,  so  most  of  it  can  be  compensated  by  a  small  axial  adjustment  and/or  tilt¬ 
ing  of  the  focusing  lens,  thus  leaving  only  a  residual  contribution  due  to  the  angular  jitter  802.  Substi¬ 
tuting  these  results  into  Eqs.  (16b,c)  one  finally  obtains 

</f“  0c)> - {X')Ay  •  V  </f0)  ( x)> ,  (19a) 


<I?Hx)>  -  y{xc„  xJAV  V„  <//0)  Gr)> 


-  V„  [xr  </J-°'(x)>  . 


One  can  obtain  expressions  for  the  relative  magnitudes  of  the  various  perturbation  terms  by 
choosing  for  Fix)  the  Gaussian  profile  defined  in  Eq.  (4).  This  yields  the  desired  target  beam 

<li0)  (x)>—  (<F>  /ir/?r)exp  (-1  x  I  2/R }),  (20) 

where  Rr  -  DT/2  —  R0  and  the  relative  perturbation  terms: 

<li"(x)>/<I±m  (3c) >  —  2  (  {xe}^  //?r>  ’  x/Rt,  (21a) 


<42’  ( x)> /  < 1 'j-0'  <x)> - [  lxcuxJAy  /Rf  +  zl2[Vli'9,JAy J  (5,4-  2 x^xJR}),  (21b) 

<I^L)  (x)>  /  <40>  (*)> - 2{  <B>2  )Ay  1  -  y  exp  (2  \x\2/2Rj) 


+  d2  (V'  <fl>  •  x! Rf, 


where  U2  -  4/-/  kRr. 


Numerical  Estimates 


If  the  alignment  and  focusing  systems  are  optimized  to  eliminate  the  average  (non-jitter)  com¬ 
ponents  of  {xc 1 4 v  and  then  the  linear  perturbations  wiil  normally  be  dominated  by  the  phase 

aberration  terms  /;  (0U  in  {xcuxc„}4y.  The  diagonal  elements  of  this  matrix,  which  predominate 
under  random  phase  conditions,  are  positive  and  will  tend  to  broaden  <  IT  (x)  > .  Their  contribution  to 
expression  (21b),  of  magnitude 


Qo  —  {\xc\2) Ay/ Rf  —  fl  9rms  / Rf  (22) 

(where  9rMS  =  ( 1 0 1 2] Vy)  indicates  that  this  perturbation  will  be  <10%  if  Rr  is  chosen  >  3 fx9 RiVf5 . 
The  magnitudes  of  the  residual  terms  in  expressions  (21a,b)  can  be  estimated  by  noting  that  the  angu¬ 
lar  jitter  89  (X'),  of  RMS  value  8 9RAts  and  transverse  scalelength  Sj,  will  have  Nj  =  Dlls}  indepen¬ 
dent  contributions  within  aperture  D2-  One  then  obtains 

Q\  —  /:  5 Orms/Rt  -Y }n  **  *y  fi  SOrmsIDi  Rr 


Qz  —  truilV  '9-<\aV\  — *  :iz&9rms/ Sj  N}n  —  C|2  89rms/D2, 


or  in  comparison  to  Q0, 


QJQo  ~  (sj/Dz)  ( Rr/fi9RMs )  (&9 r.usI 9 RM5) 


(23a) 


Qz/ Qo  ~  (z\zJ f  z)  (Rt/D2)  (Rrlfz  9rms )  (&9rMS/9rms)  (23b) 

For  typical  cases,  where  RTI /i^r.vs  ~  3-5  (corresponding  to  Q0  —  0.1—0.04), 
Rr/ Dj  —  0.01,  S//D2  <  0.05,  and  Z\z) f  <  5.  the  ratios  (23a,b)  remain  small  if  8 9RMS  <  9RMS -  The 
small  magnitude  of  (23b)  confirms  the  earlier  statement  that  the  free  propagation  path  will  have  little 
effect  as  long  as  the  coherence  zones  remain  small  in  comparison  to  aberration  scaleiengths. 

The  nonlinear  perturbations  will  normally  be  dominated  by  the  2(<5>2}i4k  term  of  expression 
(21c)  .  The  largest  contributions  to  <B>  are  likely  to  be  generated  by  optical  elements  such  as  focus¬ 
ing  lenses  or  windows  located  in  the  high  intensity  decoded  beams  beyond  the  final  amplifier  stage. 
For  example,  a  0.25  ^m  beam  of  intensity  400  MW/cm2  (  e.g..  2  J/cnr  in  a  5  ns  pulse)  will,  according 
to  expression  (9),  produce  <3>  =*  0.15  radian  in  5  cm  of  fused  quartz.  Allowing  for  an  additional 
contribution  of  0.05  radian  from  the  amplifiers,  one  finds  <3>  =  0.20  radian,  and  thus 


2[<R>;]4(,  =  8%.  The  ambient  atmosphere  within  the  long  (50-100  m)  decoding  paths  may  also  pro¬ 
duce  a  significant  contribution  to  <B>,  unless  some  care  is  taken  to  minimize  its  nonlinear  refractive 
index.  One  possibility,  which  would  also  eliminate  any  rotational  SRS,  [16]  would  be  to  overfill  these 
paths  with  one  of  the  lighter  inert  gases.  For  example,  a  calculation  of  the  third  order  nonlinear  sus¬ 
ceptibility  [17!  of  Ne  at  0.25  ^m  (using  tabulated  oscillator  strengths  [18])  gives  n-J N  =*  10-37  esu; 
hence,  Eq.  (9)  yields  <B>  =  2.8  x  10-7  p  (Amagats)  I  (MW/cm2)  L  (m),  where  p  is  the  pressure. 
Assuming  400  MW/cm2  in  a  100  meter  path  at  atmospheric  pressure,  one  obtains  <B>  =*  0.01  radian. 
The  n->J  N  of  most  other  gases  with  lower-lying  energy  levels  (including  the  heavier  inert  gases,  such  as 
Xe  [17])  can  exceed  this  value  by  at  least  an  order  of  magnitude. 

The  nonlinear  beam  steering  term  d2[V  <3>}Ay  can  arise  from  either  a  large-scale  asymmetry 
in  the  beam,  or  a  number  of  small  random  intensity  nonuniformities.  In  either  case,  its  magnitude  will 
be  approximately  d2S< B> RMS/ D2,  where  8  <B>  RMS  is  the  RMS  nonuniformity  of  <B>.  In  com¬ 
parison  to  the  2(<B>1)4y  term,  one  obtains 

d2\{V‘<B>)4y\/2  [<R>2}^k=  (d-J2D2)  &<B>rm5/{  <B>2)4y,  (24) 

which  remains  small  in  a  typical  case,  such  as  ( <B>\4y  =*  0.2  and  d-J D2  <  0.03. 

The  theory  shows  that  the  total  intensity  perturbation  due  to  both  linear  and  nonlinear  aberrations 
can  readily  be  reduced  below  20%.  Moreover,  this  perturbation  is  to  a  large  extent  controllable.  The 
spatial  averaging  process  embodied  in  coefficients  such  as  (xc  Xc)4y  and  [<B>2)/4K  tends  to  reduce 
targe  uncontrollable  shot-to-shot  variations,  while  the  x— dependent  functions  in  </f2’  (x)>  and 
</rvt>  (x)  >,  which  are  smooth  and  well-defined,  tend  only  to  broaden  the  profile  somewhat.  If  the 
coefficients  can  be  controlled  to  within  10%,  then  the  net  random  perturbation  would  be  reduced  to 
<2%. 

EXPERIMENT 

A  demonstration  of  the  above  technique  for  generation  of  smooth  output  beam  profiles  was  per¬ 
formed  with  a  small  KrF  laser  system.  Figure  2  shows  the  experimental  setup.  The  oscillator  and 
amplifier  are  U.V.  pre-ionized  discharge  lasers  similar  to  the  ones  described  in  reference  (8),  with  gain 
volumes  of  0.3  x  2  x  30  cm3  and  l  x  5  x  100  cm3,  respectively.  Nominal  pulse  duration  of  the  gain  is 
20  ns.  The  oscillator  was  operated  within  a  low-Q  resonant  cavity  formed  by  the  high  reflectivity  flat 
mirror  .V/,  and  an  uncoated  fused  silica  flat  A 42  serving  as  output  mirror.  Because  no  frequency- 
selective  elements  are  present  within  the  cavity,  the  oscillator  bandwidth  is  expected  to  —50  cm-1  [81. 


Fig.  2  —  Experimental  setup  of  the  broadband  KrF  oscillator-preamplifier  system. 
The  degree  of  spatial  incoherence  at  aperture  A  is  controlled  by  the  iris  I). 


The  oscillator  could  be  operated  to  provide  either  a  spatially-coherent  or  an  incoherent  beam  by  adjust¬ 
ing  the  aperture  of  iris  /(.  For  coherent  operation,  /i  was  reduced  to  —1  mm,  giving  an  observed  out¬ 
put  beam  divergence  close  to  diffraction- limited.  With  /|,  fully  opened,  the  oscillator  operated  in  a 
spatially-incoherent,  highly- multimode  regime  with  a  measured  beam  divergence  —2x6  mrad2. 

Output  from  the  oscillator  illuminates  a  circular  aperture  A,  which  plays  the  role  of  VDA  in  Fig. 
1.  The  optical  elements  LUL2  and  .Vf<  project  the  light  at  A  into  its  far-field  and  expand  the  image. 
They  are  equivalent  to  the  L{  of  Fig.  1.  The  beam  is  then  amplified  in  two  passes  through  the  ampli¬ 
fier.  The  size  and  shape  of  the  beam  in  the  amplifier  depend  upon  the  operating  mode  of  the  oscillator. 
In  the  coherent  regime,  the  beam  at  the  amplifier  input  looks  like  an  Airy  pattern  with  diameter  —2 
mm.  In  the  incoherent  regime,  the  beam  has  dimensions  of  0.5  x  2  cm2. 

Lenses  Lj  and  Lt  project  the  Fourier-transform  of  the  amplifier  output  beam  to  a  diffuser  0, 
which  is  equivalent  to  the  target  plane  of  Fig.  1.  The  image  of  the  beam  at  the  diffuser  is  photographed 
from  behind  by  a  simple  open  shutter  camera.  Spatial  resolution  of  this  setup  is  —0.1  mm.  This  tech¬ 
nique  is  a  good  way  to  attenuate  the  intense  laser  beam,  and  was  extensively  tested  fI9l.  Polaroid  type 
55  film  was  used  in  these  preliminary  experiments. 

Figure  3  shows  the  beam  profiles  obtained  in  the  two  regimes.  The  coherent  beam  (Fig.  3a)  col¬ 
lected  information  from  all  the  aberrations,  partial  obscurations,  amplifier  gain  striations.  and  distor¬ 
tions  of  the  optical  system.  The  incoherent  beam  (Fig.  3b)  shows  dramatic  improvement  in  its  smooth 
output  profile.  In  addition  to  achieving  a  more  uniform  illumination  of  the  aperture,  one  can  readily 
remove  most  of  the  distortions  observed  with  the  coherent  beam  from  this  smail  laser  system  by  better 
alignment  and  higher  quality  (and  undamaged)  optical  components.  However,  these  problems  are 
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much  more  difficult  to  correct  on  a  large-scale  laser  system,  and  there  the  echelon-free  ISI  will  offer  a 
significant  improvement. 


The  output  power  level  was  ~1  MW/cm2,  thus,  the  amplifier  was  already  operating  in  the  saturat¬ 
ing  regime.  When  an  earlier  version  of  the  experiment  was  performed  with  the  image  of  aperture  A 
projected  directly  into  the  amplifier,  the  saturation  caused  strong  beam  distortion  in  a  bright  ring  near 
the  edge  of  the  beam.  This  shows  that  the  amplification  of  the  beam  in  the  far  field,  as  described  in 
this  paper,  is  an  essential  requirement  for  generation  of  a  smooth  spot  on  the  target. 

Without  any  aberrations  in  the  optical  train,  one  would  expect  Fig.  3  to  show  a  perfectly  relayed 
image  of  aperture  A.  The  aberrations  result  in  a  distorted  image  in  the  case  of  a  coherent  beam.  For 
the  incoherent  beam,  however,  the  aberrations  cause  only  a  blurring  of  the  sharp  edges,  in  agreement 
with  the  theory. 

SUMMARY  AND  DISCUSSION 

Echelon-free  ISI  appears  to  be  a  promising  technique  for  using  excimer  lasers  to  generate  the 
smooth  and  controllable  target  beam  profiles  required  for  direct-drive  fusion.  In  this  scheme,  the  inten¬ 
sity  and  phase  information  needed  to  produce  the  time-averaged  profile  </r(x)>  must  be  transported 
through  the  laser  system  by  coherence  zones  small  in  comparison  to  the  apertures  and  linear  aberration 
scalelengths.  (In  other  words,  the  angular  power  spectrum  of  the  partially-coherent  light  must  remain 
broad  in  comparison  to  the  angular  perturbations  introduced  by  the  aberration.)  Under  conditions  appli¬ 
cable  to  KrF,  this  criterion  should  ensure  that  <IT(x)>  remains  relatively  insensitive  to  both  linear 
aberration  and  amplifier  saturation  within  the  entire  laser  system.  Both  theory  and  experiment  have 
shown  that  the  aberration  will  tend  to  broaden  and  smooth  <lT(x)> ,  but  will  not  introduce  any 
small-scale  structure.  The  theory  also  shows  that  the  broadening  is  to  a  large  extent  controllable 
because  it  depends  only  upon  spatial  averages  of  the  aberrated  quantities  over  the  whole  laser  aperture. 
The  estimates  presented  here  have  indicated  that  it  should  be  possible  to  control  the  functional  form  of 
</(x)>  to  within  a  few  percent,  even  with  a  large  angularly-multiplexed  system.  In  excimer  lasers, 
where  self-focusing  is  weak  (e.g.,  <B>  <  1/2  radian),  both  this  and  the  other  modified  ISI  scheme 
(suggested  in  the  fourth  paragraph  of  the  introduction)  offer  another  important  advantage  over  conven¬ 
tional  ISI  with  the  echelons  at  the  output.  They  not  only  avoid  the  problems  mentioned  in  the  intro¬ 
duction,  but  also  allow  techniques  that  could  change  the  size  and  shape  of  </r(x)>  during  the  pulse 
in  order  to  match  an  imploding  pellet. 
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